Fluctuations provide strong selection in Ostwald ripening 
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A selection problem that appears in the Lifshitz-Slyozov (LS) theory of Ostwald ripening is reex- 
amined. The problem concerns selection of a self-similar distribution function (DF) of the minority 
domains with respect to their sizes from a whole one-parameter family of solutions. A strong se- 
lection rule is found via an account of fluctuations. Fluctuations produce an infinite tail in the DF 
and drive the DF towards the "limiting solution" of LS or its analogs for other growth mechanisms. 
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Ostwald ripening (OR) Jlj] is a fascinating and generic 
process of self-organization in a physical system far of 
equilibrium. It develops in a late stage of a first-order 
phase transition, in two or three dimensions, when a 
two-phase mixture undergoes coarsening and the inter- 
facial energy decreases subject to a global conservation 
law J|||. OR continues to attract considerable atten- 
tion both in experiment Q and in theory For a 
nonlinear physicist, the problem of OR is of great inter- 
est because of a long-standing selection problem [||J^-j|] 
addressed below. 

Lifshitz and Slyozov (LS) || and Wagner || devel- 
oped a mean-field formulation of OR, valid in the limit 
of a negligibly small volume fraction of the minority do- 
mains. In this formulation, the dynamics of the distribu- 
tion function (DF) F(R,t) of the minority domain sizes 
are governed (in scaled variables) by a continuity equa- 
tion, 
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where R c (t) is the critical radius for expansion/shrinkage 
of an individual domain, while n is determined by the 
growth mechanism. (For a review of different growth 
mechanisms see Refs. pjQj ] . Most known are diffusion- 
controlled growth, n = 1, and interface-controlled 
growth, n — 0.) The dynamics are constrained by con- 
servation of the total mass (or volume) of the minority 
domains: 



R 3 F(R, t)dR = Q = const . 
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Scaling analysis of Eqs. (|l]) and (||) yields a similarity 
ansatz F(R,t) = t~^ z $ (R t' 1 / 2 ) and R c = (t/er) 1/z , 
where z = n + 2 and a — const. Upon substitution, one 
obtains a family of self-similar DFs (formally, for every 
n > —1). Each of the DFs is localized on a finite interval 
[0, u m ] of the similarity variable u — Rt~ x l z . The self- 
similar DFs can be parameterized by er, and there is a 
finite interval of allowed values of the scaled coefficient 
a. For each solution, the average domain radius and the 
critical radius grow in time like t x l 2 , while the concen- 
tration of domains decreases like t~ 3 / 2 . However, the 



coefficients in these scaling laws are <7-dependent. The 
scaling function <&(£) has a markedly different shape de- 
pending on the value of a. It should be stressed that the 
problem of OR, as described by Eqs. (|l|) and (||) is fully 
determined (of course, if one prescribes an initial condi- 
tion) . Therefore, the scaled coefficient a is an observable 
quantity. It can be determined in a direct experiment or 
simulation by measuring, at large times, the coefficient in 
the power law for the critical radius R c versus time. The 
reader is referred to Ref. || for a detailed description of 
the family of self-similar DFs for different values of n. 

An important selection problem therefore arises. It has 
a long history [^,||-^|, and its present status is as follows. 
There is only a "weak" selection within the framework 
of the "classical" model (|l|) and (0) . The "weak" selec- 
tion rule, obtained recently ||||, is the following. If the 
initial DF F(R,0) has a compact support (0, R m ) and 
is describable by a power law A(R m — R) x in the close 
vicinity of R — R m , then it is the exponent A that se- 
lects the correct self-similar asymptotic DF. The selected 
value of parameter a is 
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(n + 2){v - 1) 



(n + 2)A + n + 5 
(n + 1)A + n + 4 



(3) 



(4) 



(see Ref. || for details). The celebrated "limiting" DFs 
obtained by LS || for n = 1, and by Wagner [|| for 
n = 0, correspond to extended initial DF or, formally, 
to A — > +oo. In this case Eqs. (||) and (||) yield the 
well-known "universal" value a = 9/4 for the diffusion- 
controlled OR, n = 1, found by LS 10). On the con- 
trary, as it is clear from Eqs. (^) and (0), for initial DFs 
with compact support and finite A, different values of the 
scaled coefficient a are obtained. Finally, for those ini- 
tial DFs with compact support that cannot be described 
by a power-law asymptotics near R — R m , convergence 
to any self-similar solution is impossible . A rigorous 
mathematical proof of these results is presently available 
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It has become clear after the analyses of Refs. [p|,pUlT| 
that in order to get strong selection, one must go beyond 
the "classical" model. In the present communication I 
report on a progress in this direction. Here is an out- 
line. I will employ a mean-field cluster formulation of 
the problem and proceed to the long-time limit, when 
only large clusters and single atoms dominate. Using the 
characteristic inverse number of atoms in a cluster as a 
small parameter, I will arrive at a Fokker-Planck (FP) 
equation for the cluster size distribution function. The 
drift term of the FP equation describes growth/shrinkage 
of clusters due to an interplay between attachments and 
detachments of single atoms, and it corresponds to the 
"classical" LS-theory. The diffusion term of the FP equa- 
tion accounts for fluctuations, and it is not present in the 
LS-theory. This term becomes irrelevant at long times, 
however, it can play a very important role. Indeed, even 
if the initial DF has compact support, the diffusion term 
produces an infinite tail in the DF. As the result, the DF 
will finally approach the limiting solution of LS Q] (or its 
analogs for other growth mechanisms). 

The mean-field rate equations of the cluster model (see, 
e.g., ]l2]-[r^|) represent a natural extension of the mean- 
field continuum models, as these equations account for 
the discrete nature of atoms: 



Ni = —2K\N{ 



Ni £ K sN s 

s>2 



N 2 \ - N s 
2— + > — 

s>3 



(5) 



N a = N 1 (K„_ 1 N a _ 1 -K a N a )-^- + ^±^. (6) 

T s T s+ l 

Here N s (t) is the cluster size distribution function (s is 
the number of atoms in a cluster), the K s are the rates 
of attachment of single atoms to a cluster of size s, and 
the t s are the inverse rates of detachment of single atoms 
from a cluster of size s. Rare events of direct inter-cluster 
coalescence (coagulation) are neglected in Eq. (^); this 
requires a small volume fraction of the "cluster phase" . 
As no external source of atoms is present in Eqs. (|^) and 
(^), these equations preserve the total concentration of 
atoms: 



s=l 



(7) 



a discrete equivalent of Eq. (0) . 

For most growth mechanisms (including the growth 
processes, controlled by diffusion and by interface), the 
attachment-detachment kinetics, combined with mass 
conservation, promotes growth of larger clusters at the 
expense of smaller ones, and this process is nothing but 
OR. Therefore, if the total concentration of atoms is large 
enough, the system undergoes coarsening: the average 
cluster size grows in time, and the total number of clus- 
ters decreases. The late time asymptotics of this pro- 
cess should reproduce OR quantitatively 0,0. At the 



coarsening stage the number of clusters with small s be- 
come very small, except for the concentration of single 
atoms Ni that remains relatively large because of the on- 
going detachment processes. Therefore, one should con- 
sider the population of single atoms separately. As far 
as descriptions of clusters is concerned, one can proceed 
to the limit of s 3> 1, treat s as a continuous variable 
and use Taylor expansion in 1/s in Eqs. (||) and (||). Es- 
sentially, this derivation follows the paper of Binder [jl2| . 
However, in order to account for fluctuations, we should 
keep the second order terms in 1/s (in contrast to the 
approach of Binder who kept such terms in his descrip- 
tion of the nucleation stage, but neglected them in the 
coarsening stage). As the result, Eq. (JsJ) takes the form 
of a FP equation: 
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where 



V s (t)=K s N 1 -l/r s and D s (t) = i^Vi + 1/r. (9) 

are the drift velocity and diffusion coefficient in the s- 
space. 

A continuous version of equation for Ni follows from 
the discrete equation (||). We can assume (and check 
a posteriori) that in the late coarsening stage there is 
a (quasi-steady-state) balance between the processes of 
attachment and detachment of single atoms by large clus- 
ters: 
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while the rest of the terms of Eq. (|^) become irrelevant. 
Treating s as a continuous variable, we obtain 



Ni(t) 



Jo 



N. ds 



f™K s N s ds 



(11) 



(Again, we can assume that the number of clusters with 
small s is very small and formally shift to zero the lower 
limit of integration over the continuous variable s.) 

The same Eq. (|ri| ) can be formally obtained if we 
multiply the both sides of Eq. fl|) by s, integrate over 
s and use the conservation law Eq. (^). Then, perform- 
ing integration by parts in the two remaining terms, we 
again arrive at Eq. (|ll|). In this derivation one should 
disregarded the boundary terms produced by integration 
by parts. As will be checked later, the "upper" bound- 
ary terms sV s N s and s(d/ds)(D s N s ) at s — » oo vanish. 
In this case the third term vanishes there automatically. 
The boundary terms corresponding to the lower limit of 
integration are assumed to be negligible compared to the 
terms that we take into account. Formally, one should 
require that N s vanishes sufficiently fast at s — > 0. (Re- 
member, that the number of single atoms Ni is described 
separately.) 
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Eqs. (||)-(|ll|) (supplemented by appropriate initial and 
boundary conditions) represent a complete set of equa- 
tions for the late coarsening stage. If one neglects the dif- 
fusion term, he recovers the Lifshitz-Slyozov- Wagner de- 
scription and corresponding self-similar asymptotics and 
scalings for large times |l2],[l3| . To make this recovery ex- 
plicit, one should specify the dependences of K s and t s on 
s. Looking for scale invariance, we should assume power 
laws: K s — K\s p and t s — as q . Now, assuming a com- 
pact cluster morphology (d-dimensional spherical "drop" , 
where d is equal to 2 or 3), I demand that Eqs. (||) and 
( pT| ) (without the diffusion term) coincide (after scaling 
down the coefficients) with Eqs. (Q). This gives a direct 
correspondence between the drift velocities V s (t), enter- 
ing Eq. (§), and V(R, t), entering Eq. (|l|). Using this 
"correspondence principle" I find that the concentration 
of single atoms Ni (t) scales like the inverse critical radius 
R c (t), while the exponents p and q must be the following: 
p = (d — n — l)/d and q = (n — d + 2)/d. In particular, 
for the diffusion-controlled growth (n = 1) one obtains 
p = 1/3 and q = in three dimensions, and p — and 
q = 1/2 in two dimensions. For the interface-controlled 
growth (n — 0) one gets p = 2/3 and q = — 1/3 in three 
dimensions and p — 1/2 and q = in two dimensions 
]l5fl . Note that, returning to the original, dimensional 
version of Eq. (|l|) (see, e.g. Refs. |l(|) and demanding 
an exact coincidence with the zero-diffusion limit of the 
cluster model, we can find the coefficients K\ and a as 
well, for every growth mechanism. 

Let us return to the case of a non-zero diffusion in Eq. 
(P. As it is seen from the second formula of Eq. (|^), 
the mapping procedure, described above, completely de- 
termines the diffusion coefficient D. Simple scaling argu- 
ments show that the diffusion term becomes irrelevant at 
long times flrjfl . However, even without going into much 
detail in Eq. (||), one can see that this term produces 
(already at t > 0!) an exponentially small tail in the 
DF, even if the DF had a compact support at t = 0. In 
essence, small fluctuations transform a strictly bounded 
DF into an extended one. It was shown already by LS 
that an extended DF approaches, for long times, the lim- 
iting self-similar solution In the language of Eqs. (||) 
and (Q), one can say, therefore, that fluctuations select 

n + 2 , fn + 2\ n+1 . . 

v Q = — — and cr = —— . 12) 
n+1 \n + 1 J 

Therefore, fluctuations provide a strong selection rule in 
favor of the limiting solutions of LS (for n = 1), of Wag- 
ner (for 7i = 0) and of its counterparts for other growth 
mechanisms. 

Now we can go back and justify the disregard of the up- 
per boundary terms arising from the integration by parts 
in the derivation of Eq. (^Tj) . For a DF with an exponen- 
tially small tail at large s, the boundary terms obviously 
vanish in the scaling regime, as a power-law increase of 



V s with s is too slow to change anything. This case is 
relevant for an initially compact DF, as fluctuations pro- 
duce an exponentially small tail. For an extended initial 
DF with a power-law tail, N s (t = 0) cx s _AI , finitcness of 
the cluster concentration, 

/ N s (t)ds <oo 
Jo 

requires /i > 2. Then, using the assumed power laws for 
K s and t s and evaluating the boundary terms sV s N B and 
s(d/ds)(D s N s ) 7 we must require the following inequali- 
ties p < 1 and q > — 1. Using the values of the exponents 
p and q determined from the "correspondence principle" , 
we see that these two inequalities are satisfied for any 
n > — 1, that is for all cases of physical interest. 

One can suggest the following physical argument sup- 
porting the strong selection rule. In the absence of 
fluctuations, a bounded DF always remains bounded 
flPUTlfl . In other words, there is exactly zero probabil- 
ity to have clusters with a size larger than some finite 
time-dependent size. On the contrary, the presence of 
fluctuations leads to a small, but non-zero probability of 
the appearance of clusters with any number of atoms. 
As the dynamics of OR is very sensitive to small changes 
in the region of the largest available clusters, the pres- 
ence of the infinite tail in the DF will ultimately affect 
the whole dynamics, driving the DF towards the limiting 
DF. Of course, as fluctuations in macroscopic systems 
are extremely small, the time necessary for the DF to 
actually converge to the selected limiting solution can be 
extremely long (if one starts from a bounded DF). In 
this case I expect that, on a (quite long) intermediate 
time scale, a self-similar DF selected by A will develop, 
and only at much later times crossover to the limiting 
DF will be observed. This crossover can happen much 
earlier if coarsening in mesoscopic systems is considered, 
where the role of discrete nature of atoms increases dra- 
matically. For example, I expect this effect to be ob- 
servable in the processes of submonolayer relaxation of 
atomic clusters on surfaces, after epitaxial deposition is 
stopped. 

Let us compare the cluster approach used in this work 
with the approach of Mullins |l7]] . Mullins accounted for 
the fact that droplets with the same radius do not neces- 
sarily have the same expansion/shrinkage rates (because 
of correlations). He generalized the classical LS-model 
by replacing the deterministic growth law for a droplet, 
R = V(R,t) (where V is given by the second equation 
in (0), by the equation (R\R) = V(R,t), where (R\R) is 
the average value of R for droplets with a given R. Then 
he inserted this relationship in the continuity equation 
(|l|) . In contrast to the cluster approach, the approach of 
Mullins does not produce a diffusion term (as it does not 
account for fluctuations related to discrete character of 
particles) , therefore it does not provide a strong selection 
mechanism for OR. 
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Finally, I briefly speculate on possible additional mech- 
anisms of strong selection. An account of fluctuations 
represents only one of possible ways of going beyond the 
"classical" LS-formulation. Various finite volume frac- 
tion effects can provide alternative ways. One such alter- 
native is rare coagulation events that can be accounted 
for already in the mean-field formulation. This alter- 
native was briefly discussed by LS already in 1961 
As the result of the rare coagulation events, a DF that 
had a compact support at t = is also expected to de- 
velop a tail which will drive it towards the limiting solu- 
tion. No quantitative analysis of this scenario is presently 
available. Another possibility involves correlation effects, 
completely ignored by any mean-field description. Here 
I should mention the work of Marder 18 who studied 
screening effects in OR and arrived at a different FP 
equation for the DF. In his analysis, the diffusion term 
results from the screening effects, rather than from fluc- 
tuations, and it is proportional to the square root of the 
volume fraction. Of course, the presence of any linear 
diffusion term in the FP equation will produce a tail in 
the DF and drive the DF towards the limiting solution 
(if the diffusion term is small enough and does not in- 
terfere in the scaling regime). The comparative role of 
these possible selecting mechanisms is obviously volumc- 
fraction-dependent. Of course, non-universal transients 
and convergence rates towards the selected DF are ex- 
pected to differ significantly in the different scenarios. 

In summary, small fluctuations provide a strong selec- 
tion rule in the problem of Ostwald ripening, as they 
drive the system, at long times, towards the limiting 
Lifshitz-Slyozov solution. 

I acknowledge useful discussions with O. Biham and 
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grant from Israel Science Foundation, administered by 
the Israel Academy of Sciences and Humanities. 
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